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ACCELERATION OF SMALL ASTROPHYSICAL GRAINS DUE TO CHARGE FLUCTUATIONS 
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ABSTRACT 

We discuss a novel mechanism of dust acceleration which may dominate for particles smaller than 
0.1 /im. The acceleration is caused by their direct electrostatic interactions arising from fluctuations 
of grain charges. The energy source for the acceleration are the irreversible plasma processes occurring 
on the grain surfaces. We show that this mechanism of charge-fluctuation-induced acceleration likely 
affects the rate of grain coagulation and shattering of the population of small grains. 
Subject headings: dust, extinction - ISM: evolution - ISM 



1. INTRODUCTION 

Dust is an important constituent of interstellar 
medium (ISM) , molecular clouds and accretion disks (see 
Whittet 2003, Draine 2009). It gets involved in many key 
processes, for instance, it controls heating and cooling of 
the ISM (see Draine 2003, Tieliens 2005), reveals mag- 
netic fields through grain alignment (see Lazarian 2007 
for a review) and interferes with the attempts to measure 
properties of CMB radiation (see Lazarian & Finkbeiner 
2003, Fraisse et al. 2009). 

Small, i.e., less than ^ 10~^ cm grains, are an im- 
portant component of the interstellar dust population, 
with a notable fraction of very small grains - Polycyclic 
Aromatic Hydrocarborn (PAH) particles, which are es- 
sentially large molecule^ (see Leger & Puget 1984). In 
what follows, due to the reasons that are explained in §5, 
we do not directly address PAH particles, but our ap- 
proach to the small grain acceleration may be extended 
to this important population of grains. 

Most properties of grains, including light extinction, 
electron photoemission, and chemical activity depend 
not only on grain chemical composition, but also on their 
sizes. In astrophysical media, these sizes are affected by 
grain-grain collisions. The minimal velocities of grains 
are determined by their Brownian motion correspond- 
ing to the temperature of the ambient gas. Large-scale 
hydrodynamic motions associated with turbulence can 
make grains move faster (see Draine 1985). 

Since most astrophysical media are magnetized and 
grains are charged the hydrodynamic treatment of ac- 
celeration is frequently not adequate. A proper treat- 
ment of the grain acceleration through the interaction of 
charged grains with magnetohydrodynamic (MHD) tur- 
bulence has been developed recently (Lazarian & Yan 
2002, Yan & Lazarian 2003, Yan, Lazarian & Draine 
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^ The information about these grains is being obtained through 
the transient heating of grains by UV photons. Indeed, after ab- 
sorbing a photon a grain gets heated to rather high temperature, in- 
ducing emission of infrared photons (sec AUamandola et al. 1989). 
The PAHs were invoked by Draine & Lazarian (1998) to explain 
the anomalous emission of observed in the range of 10-100 GHz 
(see also Hoang, Draine & Lazarian 2010). 



2004, Yan 2009). This treatment makes extensive use 
of the advances of compressible MHD turbulence (see 
Cho & Lazarian 2002, 2003) and provides the mathemat- 
ical formalism of the second-order Fermi acceleration of 
charged grains interacting with MHD turbulence. 

However, the acceleration mechanisms based on the 
MHD interaction of turbulence and charged grains ex- 
hibit acceleration rates that decrease as dust particles get 
smaller. This decrease arises from the fact that the Lar- 
mor radius of charged grains becomes smaller with the 
decrease of grain mass and, correspondingly, grains have 
to interact with smaller, i.e., less powerful, turbulent fluc- 
tuations. In addition, compressible fluctuations, i.e., fast 
modes, which were identified in Yan & Lazarian (2003) 
with the most efficient acceleration, get suppressed at the 
small scales due to plasma damping, while the Alfvenic 
mode gets inefficient for acceleration at small scale due 
to anisotropy (see Yan & Lazarian 2003 for more dis- 
cussion) . The acceleration of grains with sizes less than 
~ 10^^ cm becomes rather inefficient for most media 
discussed in Yan & Lazarian (2003), Yan, Lazarian & 
Draine (2004, henceforth YLD04). 

Are there other mechanisms of dust acceleration which 
dominate for grains smaller than 10~^ cm? One can 
expect that dust-plasma interactions may be important 
for such grains. This paper presents a novel promising 
acceleration mechanism based on grain-grain Coulomb 
collisions in the presence of grain charge fluctuations. 
This mechanism utilizes intrinsic non-equilibrium nature 
of the dusty astrophysical plasmas. In addition to acting 
on the small grains, the mechanism may drive the accel- 
eration of larger grains whenever the MHD acceleration 
mechanism is suppressed (e.g., when MHD turbulence is 
damped) . 

In what follows, we introduce major timescales char- 
acterizing dynamics of charged interstellar grains in §2, 
present the mechanism of small grain acceleration as- 
sociated with charge fluctuations during their Coulomb 
collisions in §3, calculate the acceleration for grains of a 
given size distribution in §4, analyze implications of the 
effect for interstellar phases in §5, discuss the importance 
of our results in §6, and summarize them in §7. 
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2. DUST INTERACTIONS IN INTERSTELLAR GAS 

2.1. Dust damping time 

Interactions of dust with the ambient gas present the 
primary mechanism of dissipating streaming motions of 
grains. The damping rate of translational motion arising 
from the interaction with neutral gas is essentiaUy the 
inverse time for coUisions with the mass of the gas equal 
that of a grain (Purcell 1969), 
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where m„, n„, and T„ are the mass, volume density, and 
temperature of neutrals, pd is the mass density of dust 
grains and a is their radius. 

When the ionization is sufficiently high, the interac- 
tion of charged grains with a plasma becomes important 
(Draine & Salpeter 1979). The ion-grain cross section 
due to long-range Coulomb force is larger than the atom- 
grain cross section. As a result, the rate of translational 
motion damping gets modified (Draine & Salpeter 1979). 
For subsonic motions the effective damping rate is renor- 
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Here Xi is the abundance of ion i (relative to hydro- 
gen) with mass rm and temperature Ti, x = ^'^d 
f{Z) is the grain charge distribution function (see Hoang, 
Draine & Lazarian 2010). When the grain velocity Vd rel- 
ative to the gas becomes supersonic the dust interactions 
with the plasma is diminished, and the damping rate in 
this case is renormalizcd due to the gas-dynamic correc- 
tion (PurccU 1969), 
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where Cg = \pi^rjm^ is the sound speed. 

The Larmor rate for a grain is TjJ^ = \Zi^\eB jvndC^ 
where Z^ia) is the mean grain charge and B is the mag- 
netic field strength. If tl is larger than Tdamp the effect 
of magnetic field on dust dynamics is negligible. 

2.2. Grain Coulomb interactions 

Grains can interact with each other through Coulomb 
forces if the Debyc screening length is larger than the 
distance between grains. In Figure [T] we show the ra- 
tio of the Debyc radius Ad to the separation da between 
the grain of size a and grains of a larger size, calculated 
for three different idealized interstellar phases: the cold 
neutral medium (CNM) , warm neutral medium ( WNM) , 
warm ionized medium (WIM). Table [T] presents physical 
parameters for such environments [including also refiec- 
tion nebula (RN) and photodissociation region (FDR)], 
where njj is the hydrogen density, T and Tds are the 
gas temperature and dust material temperature, respec- 
tively, n(H+),n(M+) and n(Il2) are ion hydrogen den- 
sity, ion metal density and molecular hydrogen density. 
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Fig. 1. — The ratio of the Debye screening length, Ad, to the 
distance da between grains larger than a given size a. 



TABLE 1 

Idealized Environments for the Interstellar Matter 



Parameter 


CNM 


WNM 


WIM 


RN 


PDR 


nH(cm~'^) 


30 


0.4 


0.1 


10^ 


10^ 


r(K) 


100 


6000 


8000 


100 


1000 


r<i,n(K) 


20 


20 


20 


40 


80 


XH = n(H+)/nH 


0.0012 


0.1 


0.99 


0.001 


0.0001 


XM = n(M+)/nH 


0.0003 


0.0003 


0.001 


0.0002 


0.0002 


y = 2n(H2)/nH 


0. 


0. 


0. 


0.01 


0.01 



respectively. To represent grain charge we use the charg- 
ing data from Hoang et al. (2010). 

It is clear that for grains smaller than ^ 3 x 10^^ cm 
the grains do interact electrostatically with their larger 
neighbors. As a result of charge fluctuations, the dust 
particles experience forces which can accelerate them. 

We choose to depict distances to the grains of the size 
larger than the given grain size, because such interactions 
will be most important in terms of momentum trans- 
fer. Naturally, larger grains will be still interacting with 
smaller ones and therefore there should be no sharp cut- 
off for the grains to stop being accelerated by the charge- 
ffuctuation mechanism. 

3. MECHANISM OF DUST ACCELERATION 
3.1. Physics of the process 

Astrophysical grains are known to be charged (see 
Draine & Sutin 1987). The grain charges are not con- 
stant - they fluctuate around an equilibrium value Qq = 
cZq. In a plasma without external radiation Qq is nega- 
tive, because electrons are much faster than ions (in the 
presence of sufficiently intense external UV radiation it 
becomes positive, see Goree 1994). In general, one needs 
to deal with a distribution of charges (see Weingartner 
& Draine 2001, Hoang et al. 2010). 

It is possible to demonstrate that the mutual Coulomb 
collisions between dust particles with fiuctuating charges 
do not conserve the kinetic energy of the particles. As 
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the energy is being taken from the fluctuations the kinetic 
energy of the particles increase in a stochastic way (see 
de Angehs et ah 2005, Ivlev et ah 2005, henceforth 105) 
- a process similar to second-order Fermi acceleration, 
where the role of randomly moving magnetic mirrors is 
played by the randomly fluctuating electric forces exerted 
by other particles. 

The physical mechanism behind such stochastic heat- 
ing is easy to understand. Each particle creates in its 
vicinity a fluctuating electric field, so that the superposi- 
tion of these fiuctuations creates a stochastic background 
which can be stronger than the fiuctuations induced by 
the plasma itself. An individual dust particle in such 
a random environment experiences fluctuating electric 
force and hence performs a random walk in the momen- 
tum space. This naturally results in the stochastic accel- 
eration of dust particles. 

Below we elaborate on this mechanism. We first con- 
sider binary interactions (collisions) between particles 
with fluctuating charges, and demonstrate that such in- 
teractions result in the effective particle acceleration. 
Next, we employ the kinetic approach and show that 
in large particle ensembles the stochastic acceleration in 
binary collisions can result in the exponential growth of 
the particle kinetic energy (temperature) . We stress that 
we consider collisions between grains only for illustrative 
purposes, as it is one way to show that energy of grains 
do change as a result of charge fluctuations. At the same 
time, the essential element of the acceleration process 
is the fluctuation of electric force acting on a particle, 
which occurs even if grains are not directly colliding. 

3.2. Stochastic acceleration in binary collisions 

We consider dilute dust clouds and hence only focus on 
the binary interactions between grains. Binary collisions 
can be conveniently studied in terms of the center-of- 
mass and relative coordinates. Below we consider grains 
of the same mass, rridi = rnd2{= rnd), although all results 
can be straightforward generalized for arbitrary mass ra- 
tio (Ivlev et al. 2004): e.g., for mdi/md2 ^ 1 one should 
substitute the mass of light particle instead of the re- 
duced mass {^rrid). For a pair of particles with momenta 
Pi and p2, the center-of-mass and relative momenta are 

Pc ~ |(Pi + P2) and Pr = Pi — P2, respectively. The ki- 
netic energy of the pair can be expressed via pc and Pr as 
follows: pl/2md+pl/2md = pi / vid + p^J -iiud = ec+e^. 
The center-of-mass momentum and hence the energy Ec 
are conserved during the collision (here we neglect mo- 
mentum exchange due to collisions with electrons and 
ions), whereas the relative momentum is changed, 

Pc=Pc; Pr=Pr+q- 

For constant charges e-c is obviously conserved during the 
collision, but when charges vary varies as well. Thus, 
in the presence of charge fluctuations the exchange of 
the relative momentum can be divided into the elas- 
tic and inelastic parts, q = qo -t- (5q. The elastic part 
keeps the magnitude of the relative momentum constant, 
|Pr + Qol = |Pr|, so that the energy variation after the 
collision is fcr = (pr-l-qo) ■<5q/2TOd-l-(5q)^/4md. The in- 
elastic momentum exchange ^q is generally a function of 
Pc and Pr and is determined by the stochastic properties 
of charge fiuctuations. 



In order to calculate the energy variation, we suppose 
that for different particles the charge fluctuations are un- 
correlated, so that for a pair of particles with charges 
Qi{t) = Qo + 5Qi{t) and Q2{t) = Qq + 5Q2{t) we have 
{5Qi{t + T)5Q2{t)) = 0. Assuming that \5Q/Qo\ < 1, 
the inelastic momentum exchange due to the collision can 
be calculated as (5q « ^[5Qi{t) -f 5Q2{t)\^<^Qdt, where 
iy9o(r) is the potential of the equilibrium charge Qo and 
the integration is performed along the (equilibrium) col- 
lision trajectory v{t). The resulting energy variation is 
obtained by averaging over the fiuctuations, i.e.. 



1 



4TO(i 2md 



6QWipodt / SQ'Wipodt' 



(4) 

To calculate Ss-^ we also assume \SQ/Qo\ <C 1 and there- 
fore can use the following charge auto-correlation func- 
tion (Matsoukas & RusseU 1997): 
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Here ctq is the charge dispersion and is the "charging 
frequency" (which plays the role of the scale parameter 
in the charge fluctuations spectrum). The assumption 
about small charge fiuctuations simplifies significantly 
the derivation of Eq. and also allows us to treat 

the charge as a continuous variable and hence employ 
Eq. ©. 

There are several mechanisms that can result in the 
fluctuations (Yan et al. 2004, Khrapak et al. 1999). 
In this paper, for the sake of convenience we focus on 
the most elementary charging process which is due to 
discreteness of the elementary charge. Usually, this 

with an = TT^rfeleQol and 
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'.'pi 



■\/ Aiie^ni/mi is the ion plasma frequency, 

vj^i = y/k^Ti/mi is the ion thermal velocity, and z — 
e\Qo\/ ak^Ti is the dimensionless charge of dust particle 
of the radius a (for different gases, z = 1 — 3) (Tsytovich 
1997, Fortov et al. 2005). Note that other charging pro- 
cesses can also play an important role (in particular, the 
photoelectric emission) and require separate careful con- 
sideration. 

One can identify two limiting collisional regimes asso- 
ciated with the charge fiuctuations. In the limit of rapid 
fluctuations the charges vary many times in each act of 
binary collisions between dust particles. This regime cor- 
responds to t'chTint ^ 1, whcrc T{at{vd) is the character- 
istic timcscale of the dust-dust interaction during the 
collision. By employing Eq. ([5]), we readily obtain from 
Eq. (|4]) the following limiting expression for the energy 
variation: 
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dr 



dt + 0{y-;^ 



(6) 



Equation ^ shows that in this limit each collision be- 
tween dust particles results in the increase of their kinetic 
energy (105). 

The interaction timescale Tint is determined by the 
range of the fluctuating field. In contrast to the aver- 
age field, which is screened at the plasma Debye length 
Ad, the screening of the fluctuating field occurs at the 
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so-called "plasma flux" scale Aflux (see, e.g., Tsytovich et 
al 2008). The physical meaning of Aflux is rather simple: 
While the screening of the average field is characterized 
by the timescale ^ a;"/ and occurs at the length scale 
~ VTi/ojpi = Ad, the relevant timescale for charge fluc- 
tuations is ~ and hence they are screened at the 
length scale ^ Aflux = vj-Jvch- Thus, the dust interac- 
tion timescale is Tint ^ Aflux/'^^d and therefore the condi- 
tion of rapid fluctuations is reduced to VTi/vd ^ 1, i.e., 
it is generally satisfied for the subsonic dust. 

In the opposite limit of slow fluctuations the charac- 
teristic time of the charge variations is much longer than 
the collision time. This regime is completely analogous 
to "regular" Fermi acceleration - with equal probability, 
each collision results in energy gain or loss, but the col- 
lisions resulting in the energy gain are more frequent. In 
this case, after averaging over many collisions, we also 
obtain net energy growth. Analysis of this regime will 
be presented elsewhere. 

3.3. Evolution of dust kinetic temperature 

In order to understand the effect of the energy vari- 
ation in binary dust-dust collisions [given by Eq. (|6])] 
on the mean kinetic energy of the whole dust ensemble 
one should employ the kinetic approach. The kinetics is 
described in terms of the velocity distribution function 
fd{p,t). There are two principal contributions to the 
dust kinetics - one is due to the mutual dust collisions 
and another because of dust interactions with the am- 
bient gas. For simplicity, we first consider the idealized 
situation for the second contribution and only take into 
account the interactions with neutrals of temperature T„ 
(we generalize this idealized setup for more realistic en- 
vironments in §5). 

The resulting kinetic equation has the following form: 

= S^ddfd + Stdnfdi (7) 

where St^c; and Stdn denote the collision operators (inte- 
grals) describing the dust-dust and dust-gas interactions, 
respectively. Although the analysis of Eq. ([7|) can be per- 
formed for arbitrary fdip)i for the sake of convenience 
we assume that dust particles have a Maxwellian veloc- 
ity distribution /m(p)- Then the equation for evolution 
of mean kinetic energy (kinetic temperature) k^Tdit) = 
5 / (p^/™d)/Mdp is obtained by taking the second mo- 
ment of Eq. ©, 

Td^l f — (Stddfu + Stdnfu) dp. (8) 
3 J rud 

The integrals in Eq. ([8]) are calculated separately in 
Appendices |A] and |B] Equation ()Aip yields the sec- 
ond integral (contribution of dust- neutral interactions), 
whereas the first integr al d escribing the dist-dust colli- 
sions is given by Eq. (|B4p . The latter is determined 
by the kinetic coefficient Add, which is calculated in Ap- 
pendix O Eq. dCl]). By substituting this in Eq. ([Bit 
(and taking into account that the coefficient Bdd can be 
neglected, see Appcndix[B]), we derive the following equa- 
tion for the kinetic temperature of dust (105): 

Td = ^^Td - 2r-J (Td - r„) , (9) 



where ujpd — \/^7TQ^ndJrnd is the dust plasma frequency. 

Without charge fluctuations, when the collisions be- 
tween dust particles conserve the energy, the equilibrium 
temperature of dust is determined by interactions with 
the ambient gas, so that = r„. Random charge fluc- 
tuations provide an additional energy source, and if the 
coefficient of the first (source) term in the r.h.s. of Eq. 
(O exceeds the damping rate 2r^^, then the dust tem- 
perature grows exponentially with time. 

The kinetic coefficient Add (Appendix [C]) is linearly 
proportional to the maximum scattering angle, which we 
set equal to unity. Hence, we implicitly supposed that 
there are sufficiently small impact parameters p that en- 
sure scattering at large angles, X ^ 1- However, since 
the lower bound of impact parameters is limited by the 
particle radius, p ^ a, this assumption is only valid if 
the kinetic temperature is below a certain critical value 
T". From the relation x ^ Qo/pTd (Lifshitz & Pitaevskii 
1981) we readily deduce, 

rpcr QO _ \Q0\ .rp 

a e 

so that the actual value of the maximum scattering an- 
gle is equal to T"' /Td- Therefore, the exponential tem- 
perature growth described by Eq. ^ proceeds until Td 
reaches the critical value T||''(^ Ti). At larger temper- 
atures the source term in Eq. ^ tends to a constant 
value (which is obtained by replacing Td with T^'^). The 
temperature growth crosses over to linear and eventually 
gets saturated due to the friction, which determines the 
ultimate temperature of dust. 

4. SIZE-DEPENDENCE OF KINETIC TEMPERATURE FOR 
ASTROPHYSICAL GRAINS 

Let us now determine the conditions (i) when the ex- 
ponential heating described by Eq. © can set in, (ii) 
what is the magnitude of the ultimate kinetic tempera- 
ture that can be reached by astrophysical grains due to 
the stochastic heating, and (iii) what are the characteris- 
tic timescales to reach this temperature. We shall assume 
an isothermal (T„ sa w Ti) quasineutral (n^ sa Ui) 
plasma with arbitrary value of the ionization fraction 
Ui/un. Furthermore, we naturally consider dust as the 
"test species" which does not change the global energy 
balance, and therefore do not take into account the feed- 
back the dust has on the ISM plasma. 

(i) By substituting the above defined charge dispersion 

fjQ as well as frequencies ojpd, I'ch, and in Eq. ([9]) we 
obtain the heating condition, 

> —Z{2 + ZJUnO"^. 

eui 6 

Then, by using the definition of z and taking into account 
that in an isothermal hydrogen plasma z ^ 2.5 (Tsy- 
tovich 1997, Fortov et al. 2005), we reduce the condition 
for the dust heating to the following simple formula: 

a^< 0.1X1"^. (10) 

This inequality shows that the heating occurs for parti- 
cles below a certain critical size. The mechanism of the 
stochastic acceleration is effective for collisions with par- 
ticles of size > a, so that the relevant dust density to 
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be substituted in Eq. pH)) is determined by the integral 
nj^{a) = {dnj^/da)da, where dnj^/da is the appropri- 
ate size distribution. Usuahy dnj^/da decreases rather 
steeply and therefore the critical size should have a weak 
dependence on the parameters of the ISM. 

For the estimate, let us assume the MRN size distribu- 
tion for dust (Mathis, Rumpl, & Nordsieck 1977) valid 

for the size range between a few dozens of A and a few 
tenths of /zm: 

^ = -nH^MRNfl"^-', (11) 

da 

where Amrn ~ 10"^^ cm^-^ (Draine & Lee 1984). This 
yields nd(a)/nH O.SAmrno"^-^. Taking 10"^'' g, 

Hi ^ 1 cm^'^j Ti ^ 0.1 eV, and assuming also n„ « 
TT-H, we get the critical diameter of ~ 0.3 /xm, i.e., the 
heating can already be triggered in the submicron range. 
It is noteworthy that the critical size also has a weak 
dependence on the stochastic properties of the charge 
fluctuations. For instance, from Eq. © one can readily 
deduce that the scaling of the critical size on the charging 
frequency is cx; 

(ii) In order to determine the ultimate temperature 
that can be reached, we recall that when Td ex- 
ceeds the critical value T^'', the source term in Eq. 
^ is saturated at a constant value. Hence, we have 

~ 0.1(AD/a)2(nd/n„)T7. By using the defini- 
tion of z, the critical temperature can be rewritten as 
T" = z'^{Ti/eYa. Then we finally get the following esti- 
mate for the ultimate kinetic temperature of grains: 

A4n,(a) ^ ^^^^ 

Ti Un a 

Thus, the increase of the dust temperature relative to the 
temperature of gas depends very strongly on the grain 
size. Assuming again the MRN distribution for nd{a) we 
get T^/T, cx r^j^g corresponding RMS velocity of 

grains, = \J k-^T^ / rud, scales as cx a~^-^^, whereas 
the dependence on the plasma parameters is weaker, cx 

T^^'^n^ ^^"^ . For the parameters used above, a 30 nm dust 
particle can reach the kinetic temperature which is 5-6 
orders of magnitude higher than Ti. 

(iii) The heating timescale Too (i-C-, when is 
reached) is determined by the stage of linear temper- 
ature growth. This can be straightforwardly obtained 
from Eq. (O, 

°° VQo/ niirinVT,' 

The heating timescale is independent of the plasma den- 
sity and has rather weak dependence on other param- 
eters. Too fx an^^T'j This estimate yields heating 
timescale of several Myr for 30 nm grains. 

The formalism above is rather general and is applica- 
ble to any astrophysical media, including accretion disks, 
circumstellar regions and the ISM. In the next section we 
focus our attention on the ISM dust acceleration. 

5. IMPLICATIONS FOR^ ISM DUST 

The ISM dust acceleration is extremely important pro- 
cess for shattering and coagulation of dust (see Hirashita 
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Fig. 2. — RMS velocity of dust grain for the parameters of ISM 
phases given in Table 1. The grain size distribution function from 
Weingartner & Draine (2001) is adopted. Diamond symbols denote 
the transition from supersonic to subsonic motion of dust. 

& Yan 2009). We discuss both subsonic and supersonic 
motion of dust. In Eq. (|D4p (Appendix D) we rewrite 
Eq. in terms of the RMS velocity of accelerated 

grains assuming the MRN dust size distribution, Eq. 

oo 6x10^ f ( a ^-3.25 

""^^ ~ J^^^ \\mK ) VlO-6 cm J 

(13) 

where where rii is in cm ^ and the rcnormalizing factor 
a for the damping rate in the subsonic and supersonic 
regimes is given by Eq. ^ and ([3]), respectively. Note 
that in the latter case a is the function of dust velocity 
and hence Eq. (fT3|) should be resolved for . 

In Figure[2]wc compare the mechanism presented in the 
paper with the calculations of the velocities of interstellar 
grains obtained in YLD04. We see that for sufficiently 
large grains the acceleration due to MHD fluctuations 
is stronger. At the same time, the new mechanism of 
acceleration associated with charge fluctuations is more 
efficient for a < 10^^ cm. 

The efficiency of the charge-fiuctuation-induced accel- 
eration increases for smaller grains. However, the limi- 
tation of our treatment above is that the grain charging 
was considered as a continuous process [see Eq. ([5])]. 
Formally, this approach is not applicable to the situation 
when the mean grain charge Qq is so small that each 
act of losing or acquiring an elementary charge results 
in substantial discrete fluctuations. However, we believe 
that even in this case we may get right order of mag- 
nitude estimate of grain acceleration. Further research 
with Monte-Carlo simulations instead of using Fokker- 
Plank approach should test the latter conjecture. 

Comparing the Figures [T] and [5] we see the limited na- 
ture of our quantitative calculations. For grains larger 
than ~ 3 X 10~^ cm the interactions decrease due to the 
less efficient momentum deposition as grains mostly feel 
the electric field of grains smaller than their size. On 
the other hand, the charge fluctuations become large for 
grains appreciably smaller than this size and this makes 
the treatment of the problem with the Fokker-Planck 
equation not precise. Nevertheless, even taking into ac- 
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count the pilot nature of our estimates it is clear that the 
process of acceleration is potentially very important for 
small grains. 

6. DISCUSSION 

We showed that there exists a novel powerful mecha- 
nism of charge-fluctuation-induced acceleration operat- 
ing in astrophysical plasmas. This mechanism, which 
occurs to be more efficient for smaller grains, is based 
on the charge fluctuations that grains experience during 
the mutual Coulomb collisions. As the result, the kinetic 
energy of grains is not conserved in the collisions, which 
causes the stochastic heating analogous to the second- 
order Fermi acceleration. 

This acceleration mechanism is a generic plasma pro- 
cess which can operate in very different environments 
ranging from the ISM to laboratory gas discharges. It is a 
remarkable manifestation of the intrinsic non-equilibrium 
(energetic openness) of dusty plasmas (Tsytovich 1997, 
Fortov 2005): Even if electrons and ions themselves were 
in detailed equilibrium, their absorption and the subse- 
quent recombination on dust grains cannot be balanced 
by the corresponding inverse process (of the ionization 
and emission from grains) due to its apparent inefficiency 
- temperatures necessary for that should be enormously 
high. Therefore, the process of dust charging goes only 
in one direction, resulting in the imbalanced plasma flux 
on the grain surface. The dust dynamics gets coupled to 
the charging process due to charge fluctuations (which 
are inevitable due to the discreteness of the elementary 
charge). Thus, in terms of the energy balance, the mech- 
anism of charge-fluctuation-induced acceleration is based 
on the conversion of the energy flux (associated with the 
plasma flux on grains) into the kinetic energy of dust. 

One should point out that for typical ISM conditions 
the situation is very different from the equilibrium one: 
The plasma temperature is much larger that the dust ma- 
terial temperature (Tdm) due to efficient radiative cooling 
of the dust. Moreover, the ions and electrons are usually 
produced by the ISM UV and/or cosmic ray ionization. 
Grains themselves emit photoelectrons interacting with 
UV photons. All these non-equilibrium effects can alter 
charge fluctuations and therefore additionally contribute 
to our acceleration mechanism. 

In particular environments, e.g., within quiescent dark 
clouds and protoplanetary disks with suppressed turbu- 
lence, the mechanism we discussed in the paper may be 
dominant for grains substantially larger than 10~'^ cm. 
However, we notice the rapid decrease of the mechanism 
efflciency with the grain size. 

The treatment presented in this paper is limited in two 



respects. First, the Fokker-Planck approach presumes 
relatively small energy variation ([6]) occurring in each 
Coulomb collision due to charge fluctuations. If this vari- 
ation becomes large, such approach (describing the evo- 
lution of the velocity distribution function) is no longer 
applicable. On the other hand, the mean energy evolu- 
tion - which is essential for this paper - is still described 
by Eq. (|9]). Second, the direct mechanical encounters of 
grains cannot be treated within the formalism. When 
such collisions start dominating (at large Td) the ad- 
ditional dissipation associated with grain deformation, 
shattering, etc. should be taken into account. 

In future, we plan to extend the mechanism for the 
PAH population of grains, explicitly taking into account 
the substantial variation of grain charge in the process 
of individual collision with an electron or ion (remind 
that in this paper we assumed small charge fluctuations, 
which allowed us to obtain analytically tractable results). 
The fact that the efficiency of grain acceleration increases 
with the decrease of grain size makes the process very 
efficient. Extrapolating our results to PAH particles we 
may expect their collisions to be very frequent. We note 
that the additional motivation for the PAH studies comes 
from the fact that they are responsible for a component 
of CMB foreground radiation (Draine & Lazarian 1998, 
see also Hoang et al. 2010). 

Furthermore, we plan to analyze the impact of the 
charge-fluctuation-induced acceleration on the dust co- 
agulation processes occurring in the ISM and other astro- 
physical environments, e.g. protostellar accretion disks. 
It is already clear that the new mechanism can induce 
coagulation and shattering which have not been consid- 
ered in the literature before. Together with the MHD 
turbulence acceleration, our mechanism testifles that the 
astrophysical grains achieve velocities much larger than 
those arising from the Brownian motion. 

7. SUMMARY 

Our major results can be summarized as follows: 

1. The change of the grain charge in the process of 
the grain-grain Coulomb collisions results in second-order 
Fermi acceleration. 

2. The acceleration is most efficient for small, i.e., less 
than 10~^ cm grains. 

3. The process of acceleration should be accounted for 
correct description of the evolution of the smallest grains. 

AL and TH acknowledge the support of the NSF- 
funded Center for Magnetic Self-Organizatiaon (CMSO) 
and the NSF grant AST 0507164. 



APPENDIX 

A. COLLISION OPERATOR FOR DUST-GAS INTERACTIONS 

The interaction of individual dust particles with the background neutral gas is described by the Langevin equation 
(see Van Kampen 1981), 

where is given by Eq. (jlj. The random Langevin force has zero average, (L(i)) = 0, and is properly normalized, 
(L(t)L(t -|- r)) = '2t'^^ radk^Tn5{T) , to satisfy the fluctuation-dissipation theorem. 

The formalism of the Langevin equation is equivalent to the Fokker-Planck approach which describes evolution of the 
velocity distribution function (Lifshitz & Pitaevskii 1981, Van Kampen 1981). The latter approach is solely determined 
by the first and second Fokkcr-Plank coefficients, A = {6p)/5t and B — {{6p)'^)/25t (where (. . .) denotes the averaging 



7 



over many collisions occurring over the time period St). They play the role of the mobility and diffusion coefficients in 
the velocity space, respectively, and for processes that allow stable equilibrium (e.g., for dust-neutral collisions) they 
are interrelated. The resulting collision operator for dust-neutral collisions has the following differential form (Lifshitz 
& Pitaevskii 1981, Van Kampcn 1981), 

St,„/,(p)=6,„^ (jr/<^ + |^ 

where Bdn — Tdn'm-dkBTn is the second Fokkcr-Plank coefficient for dust-neutral collisions. The integration immediately 
yields 

— Stdnfudp = -6T-„ifcB (Td - T„) . (Al) 
md 

B. COLLISION OPERATOR FOR DUST-DUST INTERACTIONS 
We employ the most general form of the collision operator for binary collisions, 

Std/(p) = J [w{p',p[; P,Pi)/(p')/(p'i)- w(p,pi; p',p'i)/(p)/(pi)]dpidp'dp;. (Bl) 

This form allows for the energy non-conservation in the mutual dust collisions and therefore docs not require the 
unitarity relation to be satisfied [which reduces Eq. (jBl[) to the canonical Boltzmann form (Lifshitz & Pitaevskii 
1981)]. Here, u;(p,pi; p',Pi) is a probability function for a pair of collid ing particles with momenta p and pi 
to acquire momenta p' and p'^, respectively, after the collision. Equation (jBip counts for all possible transitions 
(p'lP'i) ~^ (P>Pi) and then is averaged over pi. Function w can be determined by solving a mechanical problem of 
the binary scattering with given interaction between the particles. 

Assuming the Maxwellian distribution, /(p) = /m(p), it is convenient to introduce the "binary" distribution function 
of dust, i^(pc,Pr) = /m(Pi)/m(P2)- Wc also define the probability functio n W (Pc,Pr; q) = w(p, pi; p',p'i) for the 
relative momentum exchange Pr — >■ Pr -I- q and after the integration of Eq. (|Bip obtain (105), 

f f \ 1 

/ — Stdd/M^p = 7 / — W/'(Pc, Pr - q; q)-F(pc, Pr - q) - W^(pc, Pr; q)^'(pc, Pr) rfpcrfprrfq- (B2) 

J rud A J nidi J 

In a homogeneous and isotropic medium, the collision probability depends only on the absolute value of the relative 
momentum, p^, and does not depend on Pc. Thus, W = W{pi-; Sq), where we introduced the absolute value for the 
inelastic momentum exchange, 6q = |pr -I- q| — Pr- For weakly inelastic collisions with Sq pr we have 

6q^2'^. (B3) 

Pr 

Furthermore, the integrand in Eq. (|B2[) can be expanded into a series over 6q. Retaining the linear and quadratic 
terms and integrating in parts, we obtain 



/ 



P^ 1 f 

Stddfudp W / {prAdd + Bdd)fMdpr, (B4) 

TUd 2md J 



where AddiPr) = J SqWdSq and BddiPr) = 5 / {Sq)^WdSq arc analogous to the Fokker-Planck coefficients for inelastic 
dust-dust collisions (105). The coefficients can be expressed in terms of the binary collision cross section S in the 
following form: 

Add ^ J Sq UdV^dY,, 

(B5) 

Bdd = \ ! {SqYridV^d'L. 

Using Eqs ^ and (|B3[) we obtain that 5q = ©(ctq). Therefore, Add = 0((Tq) and Bdd = O(crg), i.e., for weakly 
fluctuating charges the mean energy variation is determined by coefficient Add only. 

C. CALCULATION OF Add 

In order to calculate the kinetic coefficient Add, we assume the Coulomb interaction between charged grains (the main 
contribution to Add comes from the interaction at shorter distances, see 105). By substituting dipa/dr ~ —Q^/r^ in 
Eq. ([6]) and integrating along the trajectory r{t) = p'^ + (vi-t)^, we obtain the energy variation Se^ ^ (o-qQq/ p'^i'chPr) 
(here p is the impact parameter of colliding particle s, fo r simplicity we consider the scattering at small angles) an d th e 
corresponding momentum variation Sq [using Eq. (|B3p]. Finally, we substitute the result in the first equation (jBSp . 
where we use the relation between the scattering angle and the impact parameter in the Coulomb limit, x Qo/P^n 
and employ the Coulomb differential cross section, ~ {Qo/£iX^)dx (Landau & Lifshitz 1976). By integrating over 
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< X < 1 wc derive, 



Add^ ^2 ■ (CI) 



D. NUMERICAL ESTIMATE FOR RMS GRAIN VELOCITY 
The ultimate kinetic temperature of grains, T^, is given by Eq. It is determined by the plasma Debye length, 

where is in cm~'^. Let us assume the MRN size distribution for dust, Eq. (HI]). Taking into account that the 
relevant dust density is nd{a) ~ {drid/ da)da, assuming that r7,„ ^ tt-h, and using Eq. (jDl[) we get 



Ti V ani yVlOOKy VlO-« cm 



(D2) 



where the renormalizing factor a for the subsonic and supersonic regimes is given by Eq. ([2]) and ([3]), respectively. 
The corresponding mean squared velocity of grains is 

{v¥J' = = TT X 10^ i-^y (cm/s)^ (D3) 

nid VIO cm/ 

where the dust mass rud — (47r/3)pa'^ w 10~^'^(a/10~* cm)'^ is used. By substituting (jD2[) into (|D3|) we finally get 

6 X 10^ / Tj X^^V a \-3.25 
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